A0-A127  545  COMPUTER-SIMULATION  STUDIES  OF  GRAIN  BOUNDARV  SVSTEMS  1/1 
AND  MOLECULAR  SOLIDS<U>  MASSACHUSETTS  INST  OF  TECH 
CAMBRIDGE  S  VIP  2B  APR  83  ARO-15199.  7-PH 
UNCLASSIFIED  DAAG29-78-C-8806  F/G  9/2  NL 


HE  COPY 


Unciass i f 


SECURITY  CLASSIFICATION  of  this  page  ("ho*  Dmtm  Entorod) 


REPORT  DOCUMENTATION  PAGE 


REPORT  NUMBER 


4.  TITLE  fand  Subtill*) 


Computer-simulation  Studies  of  Grain 
Boundary  Systems  and  Molecular  Solids 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


S.  TYPE  OF  REPORT  ft  PERIOO  COVERED 

Final  Report 

Jan.  1,  1978  -  Jan.  31,  1983 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHOR!*) 


ft.  CONTRACT  OR  GRANT  NUMBER!*) 


Sidney  Yip 


DAAG-29-78-C-0006 


*.  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

U.  S.  Army  Research  Office 
Post  Office  Box  12211 
Research  Triangle  Park,  NC  2770 


.  MONITORING  AGENCY  NAME  ft  AOORESS!!!  dlllerent  . 


IZ.  REPORT  DATE 

April  20,  1983 


13.  NUMBER  OF  PAGES 


Controlling  Oltleo) 


14-  DISTRIBUTION  STATEMENT  (ot  tbit  Report) 


Approved  for  public  release;  distribution  unlimited. 


-  .  '  1  r 


[17.  DISTRIBUTION  STATEMENT  (ot  tho  mbmtrmct  on torod  In  Block  20,  It  dlttoront  from  Roport) 


is.  supplementary  notes 


The  view,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the 
author(s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation 

IS.  KEY  WORDS  CConlinu*  an  teretee  aid*  It  noceeeory  and  Identity  by  block  number) 

computer  simulation,  materials  properties  and  behavior,  grain  boundary 
diffusion,  grain  boundary  phase  transformation 


120.  ABSTRACT  fCrmUmm  cm  rororoo  aM*  ft  i 


See  reverse  side 


I  Idontlly  by  block  number) 


83  05  03  009 


no  romt 

W  |  JAM  71 


EDITION  OF  I  NOV  SS  IS  OBSOLETE 


UNCLASSIFIED _ 

1FCIIRITY  Cl  ASSlFlCATlOW  OF  THIS  PAGE  !■ 


t  Dote  Emoted) 


SfCUjjlTV  CLASSIFICATION  OF  THIS  PAGCfWlan  Data  EntatadJ 


Abstract 


Structural,  thermodynamic  and  mechanical  properties  of  bicrystals 
have  been  studied  using  molecular  dynamics  and  Monte  Carlo  simula¬ 
tion.  A  grain  boundary  melting  transition  in  a  two-dimensional 
system  has  been  observed.  Atomic  diffusion  along  grain  boundaries 
was  investigated  by  following  in  detail  the  migration  of  a  vacancy 
initially  introduced  into  the  bicrystal.  The  activation  energy  for 
vacancy  migration  in  iron  has  been  determined,  and  thermal  activa¬ 
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Structural  phase  transformation  in  stressed  solids  has  been  studied 
using  an  extended  Monte  Carlo  method.  Stress-strain  relations 
obtained  for  iron  showed  the  existence  of  bcc  < — >  fee  transition 
under  tension  or  compression. 

Dynamical  properties  of  a  one-dimensional  nonlinear  lattice  (the 
U^-lattice)  have  been  analyzed  in  terms  of  kink  propagation  and 
damping  and  kink  diffusion.  Both  perturbation  theory  and  molecular 
dynamic  calculations  were  carried  out. 
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I. .  INTRODUCTION 

This  report  summarizes  the  work  performed  during  the  contract 
period,  January  1,  1978  through  January  31,  1983.  The  total  funding 
was  $222,485.  The  personnel  supported  were  one  principal  investigator 
and  six  graduate  students  who  received  five  Ph.D.  and  one  M.S.  as 
terminal  degrees. 


II.  DESCRIPTION  OF  PROBLEMS  STUDIED 

The  basic  objective  of  the  project  was  to  study  by  means  of 
the  discrete  particle  simulation  techniques  of  molecular  dynamics 
and  Monte  Carlo  several  problems  in  materials  properties  and  behavior 
of  solids.  These  problems  involved  highly  nonlinear,  nonhomogeneous , 
or  nonequilibrium  processes  which  were  difficult  to  analyze  analyti¬ 
cally  or  to  measure  experimentally.  In  each  case,  appropriate  com¬ 
puter  simulation  programs  were  developed  and  results  obtained  which 
not  only  established  the  validity  of  the  simulation  modeling  approach, 
but  also  provided  detailed  information  and  insight  at  the  molecular 
level  of  materials  behavior. 

Structural  Stability  of  High-angle  Grain  Boundary  at  Elevated  Tempera¬ 
tures  and  Point  Defect  Motions  in  Self-diffusion 

Grain  boundary  structures  are  often  determined  by  static  cal¬ 
culation  in  which  the  potential  energy  is  minimized.  We  have  studied 
by  molecular  dynamics  the  structural  stability  of  such  a  configuration 
in  the  case  of  a  E-5  (9-36. 87°)  [001]  tilt  boundary  in  bcc  iron  when 
the  system  is  heated  up  to  about  two-thirds  of  melting  temperature 
and  beyond.  The  thermal  motion  of  a  vacancy,  initially  introduced 
into  a  statically  relaxed  position,  was  observed  in  sufficient  detail 
and  duration  to  determine  the  activation  energy  for  migration  and  the 
various  jump  frequencies  to  sites  in  the  grain  boundary  core.  The 
simulation  also  provided  data  on  thermal  activation  of  Frenkel  pairs 
in  the  core  and  the  mobility  of  the  interstitials  thus  produced.  In 
another  study,  simulation  was  carried  out  on  a  two-dimensional  bi¬ 
crystal,  Z-7  (9-38. 218),  to  investigate  the  existence  of  a  grain 
boundary  phase  transition.  In  addition,  grain  boundary-point  defect 
interactions  were  studied  by  Monte  Carlo  simulation. 

Crack  Tip  Processes  and  Ductile/Brittle  Behavior  in  Mode  I  Fracture 

In  atomistic  study  of  crack  propagation,  it  is  essential  to 
treat  properly  the  boundary  between  the  interior  (discrete)  region 
containing  the  crack  tip  and  the  exterior  (elastic  continuum)  region 


where  the  stress  is  being  applied.  Existing  methods  are  not  satis¬ 
factory  because  the  boundary  is  not  sufficiently  flexible  to  allow 
nucleation  of  dislocations  at  the  tip  and  subsequent  propagation 
away  from  the  crack.  A  new  method  was  developed  which  replaces  the 
exterior  region  by  appropriate  stresses  on  the  boundary  and  which 
allows  the  boundary  atoms  to  move  freely  otherwise.  This  method  was 
incorporated  in  a  molecular  dynamics  investigation  of  crack  extension 
(mode  I  fracture)  in  bcc  iron  and  fee  copper  by  determining  the  cri¬ 
tical  stress  intensity  factor  and  detailed  stress  and  strain  distri¬ 
butions  surrounding  the  crack. 

Thermal  Ignition  of  an  Exothermically  Reacting  Fluid 

A  basic  problem  in  the  theory  of  thermal  explosion  is  the  de¬ 
termination  of  the  critical  condition  under  which  the  temperature 
of  a  piece  of  material  undergoing  self  heating  rises  precipitously. 
Using  a  molecular  dynamics  model  of  two-dimensional  fluids  of  hard 
sphere  particles  capable  of  collision- induced  exothermic  reactions, 
we  have  investigated  the  critical  properties  of  slabs  of  varying 
sizes  and  densities,  results  which  can  be  compared  with  predictions 
of  classical  continuum  theory.  Other  results  include  the  induction 
period  and  temperature  profiles. 

Propagation  and  Diffusive  Motions  of  Kink  Structure  in  a  One-dimen¬ 
sional  Nonlinear  Lattice 

A  linear  chain  of  atoms  coupled  to  nearest  neighbors  by  har¬ 
monic  forces  but  each  is  also  subjected  to  a  double-well  potential 
(the  U^- lattice  model)  can  exhibit  highly  nonlinear  dynamical  behavior. 
In  particular,  there  exist  atom  configurations  where  atoms  oscillating 
in  the  left  side  of  their  double  wells  are  separated  by  a  "domain 
wall"  or  "kink"  from  those  oscillating  in  the  right  side  of  the  poten¬ 
tial  wells.  We  have  studied  the  propagation  and  damping  of  such  a 
kink  in  various  velocity  regimes  using  molecular  dynamics  simulation 
and  perturbation  theory  analysis.  Simulation  data  on  kink  diffusion 
in  a  thermal  lattice  have  been  obtained  at  various  temperatures. 


Structural  and  Mechanical  Properties  of  Stressed  Crystals  and 
Bicrystals 


We  have  extended  the  conventional  Monte  Carlo  method  to  the 
simulation  of  stressed  solids.  Structural  phase  transitions  which 
occurred  in  cubic  crystals  under  uniaxial  tensions  or  compressions 
have  been  observed  through  a  detailed  mapping  of  stress-strain 
curves  for  bcc  iron  and  fee  Ni.  Similar,  but  less  extensive,  results 
were  obtained  for  bicrystals  of  a  noble-gas  element. 


III.  SUMMARY  OF  MOST  IMPORTANT  RESULTS 


In  each  of  the  areas  of  investigation  described  above  we  have 
demonstrated  that  atomistic  simulation  can  indeed  provide  unique  and 
valuable  information.  The  most  significant  results  that  emerged  are 
the  following: 

1.  Observation  of  vacancy  jumps  in  a  bicrystal  of  bcc  iron  and 
demonstration  that  vacancy  exchange  is  the  dominant  mechanism 
for  fast  diffusion  along  grain  boundaries.  Related  results  are 
the  determination  of  activation  energy  for  vacancy  migration  in 
bcc  iron  and  the  observation  of  thermal  activation  of  boundary 
interstitials  which  were  found  to  be  largely  immobile  [1,2].* 

2.  Observation  of  grain  boundary  phase  transformation  in  a  two- 
dimensional  bicrystal  by  molecular  dynamics  simulation  at  con¬ 
stant  pressure  [8]. 

3.  Development  of  a  new  boundary  condition  for  modeling  a  system 
of  discrete  atoms  surrounded  by  an  appropriately  stressed 
elastic  medium,  and  demonstration  that  crack  extension  occurs 
in  bcc  iron  without  dislocation  emission  and  in  fee  copper 
with  dislocation  emission  [3]. 

4.  Demonstrated  ability  to  simulate  thermal  ignition  processes  in 
reacting  fluids  and  quantitative  test  of  classical  continuum 
theory  of  thermal  explosion.  Feasibility  of  studying  effects 
of  fluctuations  in  the  vicinity  of  criticality  is  established 
[4]. 


Numbers  in  brackets  refer  to  published  work  in  the  List  of 
Publications. 


6.  Determination  of  grain  boundary  interactions  with  vacancies 
and  impurities. 

7.  Development  of  a  new  dynamical  theory  of  kink  propagation  and 
damping  in  one-dimensional  nonlinear  lattices  and  study  of 
kink  diffusion  by  molecular  dynamics  simulation. 
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